Exotic order in simple models of bosonic systems 
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We show that simple Bose Hubbard models with unfrustrated hopping and short range two-body 
repulsive interactions can support stable fractionalized phases in two and higher dimensions, and in 
zero magnetic field. The simplicity of the constructed models advances the possibility of a controlled 
experimental realization and novel applications of such unconventional states. 



Recent theoretical developments [jj], |^, ||, [Ij have shown 
that two or three dimensional strongly correlated systems 
in zero magnetic field, could display quantum phases with 
fractional quantum numbers. This theoretical progress, 
inspired mostly by the search for a theory of the high- 
temperature superconductors ||, is likely to play an im- 
portant role in our eventual understanding of the mys- 
terious properties of several strongly interacting elec- 
tronic systems. However, to date, no such experimental 
system has been unambiguously shown to display frac- 
tional quantum numbers. Further impetus for the search 
for experimental realizations of fractionalization comes 
from the possibility of using such states to construct 
qubits|| |7]. The topological structure inherent in these 
states naturally protects the system from decoherence. 

The primary goal of this paper is the identification 
and possible design of specific condensed matter sys- 
tems which display the phenomenon of fractional quan- 
tum numbers. To that end, we study particularly sim- 
ple models of bosons with unfrustrated hopping and 
short ranged two-body repulsive interactions on a two- 
dimensional (2D) square lattice. We show that in partic- 
ular parameter ranges, a fractionalized insulating phase 
exists where there are excitations whose charge is one 
half that of the underlying bosons. Superfiuid or more 
conventional insulating phases result in other parame- 
ter ranges. The simplicity of our models opens up the 
possibility that they can be realized in arrays of quan- 
tum Josephson junctions, or possibly in ultra-cold atomic 
gases. This would provide a definite experimental realiza- 
tion of a fractionalized phase which could then possibly 
be exploited to construct topologically protected qubits. 

The fractionalized phase appears in a region of inter- 
mediate correlations where neither the boson kinetic en- 
ergy nor repulsive potential energy completely dominates 
over the other. This lends support to the general notion 
that fractionalization is to be looked for in a many-body 
system at intermediate correlations. For example, in the 
interacting electron system, fractionalization possibly oc- 
curs at intermediate values of density somewhere between 
the extreme low density Wigner crystal and the high den- 
sity Fermi liquid regimes. Similarly, electronic Mott in- 
sulators that are close to the metal-insulator transition 
may be good candidates for fractionalization. 

The generalization of our models to three dimensions 



(3D) is of some interest. The 3D version of our boson 
Hubbard model has in fact two distinct fractionalized 
insulating phases: First, there is a fractionalized phase 
similar to the one in 2D, with the distinct excitations 
being a charge-1/2 chargon and a Zi vortex (vison) line. 
The topological order in this phase is stable up to a fi- 
nite non-zero temperature. Experimental realization of 
this phase may therefore be of interest for the quantum 
computing application as a way of controlling errors due 
to non-zero temperature. Another distinct fractionalized 
insulator also appears in 3D. In this phase, the excita- 
tions are a gapped charge-1/2 chargon, a gapless linear 
dispersing "photon" , and a gapped topological point de- 
fect (the "monopole"). Wen[^| has recently pointed out 
that stable mean field theories may be constructed for 
quantum phases where a massless U(l) gauge boson (a 
photon) emerges in the low energy description. Our re- 
sults provide an explicit and concrete model for such a 
phase. 

Fractionalization of bosons in two dimensions: 

Consider bosons moving on the lattice shown in Fig. [l]. 
A physical realization may be a Josephson junction array 
with superconducting islands arranged on the sites of the 
"bond-centered" square lattice and Josephson-coupled 
with each other as indicated by the links. We also stipu- 
late repulsive interactions between the bosons ( "charging 
energy") that favors charge neutrality not only on indi- 
vidual islands but also on the shaded clusters (note that 
neighboring clusters share one site)^). The correspond- 
ing Hamiltonian is 
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Here, b\ = e l6r represent bosons (Cooper pairs) residing 
on the corner sites of the lattice, and ifr , 



rep- 



resent bosons on the bond-centered sites (identified by 
the bond end-points); n b r , , are the corresponding bo- 



son numbers, [9 r ,n r 



and similarly for the ?/>-bosons. 



Throughout, we work with a number-phase (quantum 
rotor) representation of the bosons, as is particularly ap- 
propriate in the Josephson junction array realization. 

The wi term is a boson hopping (Josephson coupling) 
between the corner and the bond-centered sites, and 
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FIG. 1: Josephson junction array on a 2D bond-centered 
square lattice modeled by the Hamiltonian Eq. (|l]). Each 
shaded area indicates schematically cluster charging energy 
UN?. 



r' e r sums over all such bonds emanating from r. 
The W2 term is a boson hopping between the neighbor- 
ing bond-centered sites as indicated with dashed lines in 
Fig. |} The Ub and terms represent on-site boson re- 
pulsion, while the U term is the cluster charging energy 
that favors charge neutrality in each cluster. The opera- 
tor N r associated with each cluster is defined through 

JV r = 2n* + (2) 

r'Gr 

The total boson number of the system is iV" to t = h J2r 
Both the 6-bosons and the V'-bosons are assigned charge 
qb. The model has only a global U(l) charge conservation 
symmetry. 

For large W\,W2 3> Ub,u^,U the system is a superfluid. 
In the opposite limit, Ub,u^,U ^> 11)1,11)2, the system is 
a conventional Mott insulator with charge quantized in 
units of qb- We argue below that when the charging ener- 
gies U and Ub,u^, are varied separately, there is an inter- 
mediate regime !7 > wi,W2 > y/ubU, yJu^U, in which 
the system is a stable fractionalized insulator with charge 
^ excitations and charge visons above a ground state 
with no conventional broken symmetries. A schematic 
phase diagram of our model is shown in Fig. |^. 

The analysis in the limit of large cluster interaction 
U Wi, W2, uj,, u.0 is similar to that in the large U limit 
of the electronic Hubbard model at half-filling. If the 
other terms are all zero, there is a degenerate manifold 
of ground states specified by the requirement N r = 
for each r. This ground state sector is separated by a 
large charge gap U from the nearest sectors. Including 
the wi,W2,Ub, u^, terms lifts the degeneracy in each such 
zeroth-order sector, and this is best described by deriv- 
ing the corresponding effective Hamiltonians for small 
perturbing couplings. 

Consider the ground state sector N r = for all r. An 
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FIG. 2: Schematic phase diagram of the boson Hubbard 
model Eq. (Q) for a particular cut w = wi ~ W2 and uj, ~ 
through the parameter space. 
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FIG. 3: Diagram of the large U insulating phases of the boson 
model Eq. Qlt) in two dimensions. The effective Hamiltonian 
Eq. (j^) is equivalent to the (2 + 1)D compact QED gauge 
theory coupled to a charge 2 scalar. The Mott insulator is 
conventional or fractionalized depending on whether the ef- 
fective gauge theory is confined or deconfined. 

elementary calculation gives 

H eS = Hub,^ - -'"bond ^ \i^l rl ) 2 b r br' + k.C. 
(rr') 

-i^ring ^2 (^12^23^34^41 + h.C.J , (3) 

□ 

where H UitU , stands for the on-site repulsion terms as in 
Eq. (g), Jbo„d = w\/U, and K ling = 2w%/U. 

A simple change of variables shows that to- 
gether with the constraint N r = can be regarded as 
the well-studied [jlO| (2 + 1)D compact [7(1) gauge the- 
ory coupled to a charge 2 scalar field. In (2 + 1)D, 
there are two distinct phases shown in Fig. ^|. For 

-'bond; -^"ring 

< Ub,u^, the gauge theory is "confined", 
and all excitations carrying non-zero "gauge charge" are 
confined. Zero gauge charge excitations carrying physi- 
cal charge quantized in units of qb of course exist with a 
gap of order 2U. This is the conventional Mott insulator 
of our boson model. 
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In the opposite regime, Jbond, ^ring ^ Ub,u^,, the 
gauge theory is in the "deconfined Higgs" phase. Ob- 
jects with N r — 1 at some site, i.e. physical charge ^ 
(chargon), have gauge charge 1, are not confined and can 
propagate above a finite gap of order U. There is also 
a stable gapped Zi vortex excitation (vison). The de- 
confined phase has a topological order f0], Q |: e.g., the 
ground state is two-fold degenerate on a cylinder, ob- 
tained by threading no or one vison through the hole of 
the cylinder. 

The details of the chargon motion are determined by 
the effective Hamiltonians that obtain in the charged sec- 
tors. Straightforward calculation shows that the presence 
of the chargon induces a weakening of the background 
on the bonds and plaquettes that are connected to it: 
J bond = i^bond, i^ing = faring- Of course, chargon 
confmcmcnt/dcconfincment is controlled entirely by the 
bulk Jbond) -faring vs Wf,, Uty terms that obtain far away 
from the chargon location and is as expected by looking 
at H$ only. 

Ground state wavefunctions and topological or- 
der: A good caricature for the ground state wavefunction 
of , Eq. (||), is obtained by "Gutzwiller" -projecting 
a superfluid state into the sector N r = 

|$) = Po \e r = <j> rr , = 0) = £'|K,nt», (4) 

where the last form is written in the boson number basis 
and the primed sum is over all configurations such that 
N r = at every site r. This is the exact ground state 
wavefunction when Ub = — but is not normalizable. 
A normalizable wavefunction is obtained by introducing 
a cutoff for large occupation numbers at each site as is 
appropriate for non-zero it;,, u^. Below we leave any such 
cutoff procedure implicit. 

A topologically distinct ground state on a cylinder is 
obtained by Gutzwiller-projecting a superfluid state with 
one vortex threading the cylinder |13) : 



1$, 



(5) 



where N^ ol is a sum of nf r , in a given columnar "cut" 
of x-directed links (assuming the cylinder is defined by 
periodic boundary conditions along x). Due to the con- 
straints, the parity of N^ ol is the same for all columns so 
that the location of the cut is arbitrary. The projected 
vortex state describes one vison threading the hole of the 
cylinder. 

The presence of topological order is established by 



$„)/($|$) goes 
13], and that all 



noticing that the normalized overlap (<f> 
to zero as 0(e~ cLv ) with the system size 
local physical operators are the same in the two states 
since the column defining N^, can be deformed away 
from any such operator. Thus, the states with no or one 



vison are indeed orthogonal to each other and degenerate 
in the thermodynamic limit. 

The boson Hamiltonian Eq. ([!]) is unfrustrated, in the 
sense that the hopping amplitudes are all positive. It is 
well-known then that the ground state wavefunction is 
unique and has positive amplitudes in the boson number 
basis. This does not contradict the topological order in 
the fractionalized state. The column parity (— 1) ooi is 
conserved by the Hamiltonian H^g , and the theorem ap- 
plies to H^g only separately in the even and odd sectors 
on the cylinder. By taking the combinations |<£>) ± |$„) 
we indeed obtain positive wavefunctions that reside com- 
pletely in the even or odd sectors. As far as the bare 
boson Hamiltonian Eq. (]]]) is concerned, it is more ap- 
propriate to speak of the states with no or one vison. 

Variations in 2D: A simple variation of the model 
considerably enhances the region of stability of the frac- 
tionalized phase. Consider "half-filling" for the site 
bosons described by the modified Hubbard repulsion 
terms 



H Ub: u — u b 



E 



nl- 



(6) 



All other terms are unchanged. In the large U 
limit, the corresponding compact QED theory now has 
static charges ±1 placed on the A and B sublattices 
respectively j|, and is at half- filling for the gauge charge 
2 matter field. The fractionalized insulator now occupies 
a larger area and extends all the way to infinitesimally 
small Jbond/itfe for infinitely large K V i ng /u^p as shown in 
Fig. ||. It is also more stable for large Jbond /ub due to 
frustration coming from the Berry phase terms in the 
corresponding Ising gauge theory pl|. 

In the model Eq. (|l|), chargons are bosonic excitations. 
However, as noted above, the background couplings are 
weakened in the vicinity of a chargon. It is then plau- 
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FIG. 4: Large U insulating phases of the 2D boson model at 
half-filling, Eq. (|^). The region of stability of the fractional- 
ized phase is enhanced compared with Fig. ^ at integer filling. 
There may be several (?) non-fractionalized insulating states 
with broken translational invariance; here, we focus on the 
fractionalized state only. 
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sible that the corresponding K' tin& < in some model, 
in which case it is energetically favorable for a vison 
to bind to the chargon thus forming a fermionic exci- 
tation. This possibility is indeed realized when we mod- 
ify our model slightly by allowing some frustration in 
hopping Jl4[. Thus, we obtain an explicit boson model 
that has an insulating phase with fermionic excitations 
carrying fractional charge. If this unusual insulator is 
doped, we would obtain a Fermi liquid, i.e., a metallic 
phase in a boson model, a true Bose metal. 

Three dimensions: Consider now the 3D version of 
our boson model on a bond-centered cubic lattice. Pro- 
ceeding as before, the large U Mott insulating states 
are described by the effective Hamiltonian H^g , Eq. (||), 
which is equivalent to the (3 + 1)D version of the com- 
pact U(l) gauge theory coupled to a charge 2 scalar. The 
phase diagram is shown in Fig. g, and now has three dis- 
tinct phases. The "confined" phase is the conventional 
Mott insulator. The "deconfined Higgs" phase is simi- 
lar to the fractionalized phase in two dimensions. The 
distinct excitations here are gapped charge % chargon 
and neutral Z2 vortex loop (vison) . The vison excitation 
energy is proportional to the loop length, and the loops 
do not proliferate up to a finite temperature. Thus, in 
this phase in 3D, the topological order is stable for small 
finite temperature (unlike 2D where a finite density of 
thermally excited point visons destroys the topological 
order at any non-zero temperature). 

Finally, the "Coulomb" phase is also fractionalized. In 
this phase, the low-energy theory in the ground state sec- 
tor is that of the pure gauge (3 + 1)D compact QED in 
its Coulomb phase, and has a gapless linearly dispersing 
gauge boson (photon) and a gapped topological point de- 
fect (monopole) as its distinct excitations. The charged 
sector has charge % excitations above the gap U, but 
thpsp nnw interact aria an pm.prgpn.t. kmgzrange Ckmlomb 
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FIG. 5: Same as in Fig. ^ but in three dimensions. There is an 
additional fractionalized phase, the "Coulomb" phase, with 
the distinct excitations being a gapped charge % chargon, 
gapless photon, and a gapped monopole. 



system even when the system is cycled from the super- 
fluid to the fractionalized insulator and back. For this 
to work, one needs to go to temperatures well below the 
vison gap. Such flux trapping corresponds directly to the 
ability of the qubit to retain its state in the experimental 
environment. In this context, we want to note again that 
in 3D, unlike 2D, the topological order is not destroyed 
by small finite temperature. 
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